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CONSISTENT ESTIMATES OF DEFORMED ISOTROPIC 
GAUSSIAN RANDOM FIELDS ON THE PLANE 

By Ethan Anderes^ and Sourav Chatterjee^ 

University of California at Davis and University of California at Berkeley 

This paper proves fixed domain asymptotic results for estimating 
a smootfi invertible transformation / : ^ when observing the 
deformed random field Z o / on a dense grid in a bounded, simply 
connected domain SI, where Z is assumed to be an isotropic Gaus- 
sian random field on R^. The estimate / is constructed on a simply 
connected domain U, such that U C O. and is defined using kernel 
smoothed quadratic variations, Bergman projections and results from 
quasiconformal theory. We show, under mild assumptions on the ran- 
dom field Z and the deformation /, that / — > Ref + c uniformly on 
compact subsets of U with probability one as the grid spacing goes to 
zero, where Re is an unidentifiable rotation and c is an unidentifiable 
translation. 

1. Introduction. The use of deformations to model nonstationary pro- 
cesses was first introduced to the spatial statistics literature by Sampson 
and Guttorp [37]. Their work, as well as that of subsequent authors (see, 
e.g., [13, 23, 34] and [38]) consider estimating the deformation / when ob- 
serving a deformed random field Z o / at sparse observation locations with 
independent replicates of the random field. 

Two recent papers, [10] and [5], study the different problem of estimating 
a deformation / from dense observations of a single realization of a deformed 
isotropic random field Zo / in two dimensions. These deformed isotropic ran- 
dom fields provide a flexible semi-parametric model of nonstationarity for 
random fields. In addition, this observation scenario is becoming increas- 
ingly important with the abundance of high resolution digital imagery and 
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remote sensing. One of the more recent motivations for the deformation 
model under the one-reahzation observation scenario is gravitational lensing 
of the cosmic microwave background (CMB). The gravitational effect from 
intervening matter distort the CMB images to produce deformed random 
field observations. In the hope of improving estimates of cosmological pa- 
rameters and the mass distribution in the universe (including dark matter) 
there is considerable interest in detecting and measuring the lensing of the 
CMB (see, e.g., [22] and [40]). 

In this paper, we establish a strong consistency result for the estimation 
of the deformation / when observing Z o / on a dense grid in a bounded, 
simply connected domain in M^, as the grid spacing goes to zero. We first 
construct estimates of the complex dilatation and log-scale of the map / (see 
Section 4), which converge uniformly on compact subsets of the observation 
region with probability one. Then, we construct a deformation estimate / 
on a subset of the observation region that converges uniformly on compact 
subsets with probability one. We show this result under mild assumptions 
on the map / and the two dimensional isotropic random field Z. 

Most attempts at recovering the deformation / from a single realization 
of Z o / rely on estimating local properties of /, usually related to the Ja- 
cobian of /, from the local behavior of the random field Z o f. Intuitively, 
the random field Zo / is locally stretched and sheared by /, as determined 
by the Jacobian. One can clearly see the visual consequences of this shear, 
as seen in the left plot from Figure 1. When the random field Z is isotropic, 
the identification of all four parameters of the Jacobian becomes difficult 
from the local behavior of Z o f. In particular, by decomposing the Jaco- 
bian matrix as U AV'^ (using singular value decomposition so that U, V are 
orthogonal matrices, and A is diagonal) the rotation matrix U becomes par- 
ticularly hard, if not impossible, to estimate when observing Z o / in a small 
neighborhood. An important object for us is the complex dilatation and 
log-scale of /, determined by the Jacobian (and defined in Section 4), which 
are invariant under left multiplication of rotation matrices to the Jacobian. 
It is this invariance that allows us to estimate these parameters under the 
isotropy assumption for Z. 

Guyon and Perrin [20] tackle the problem of developing consistent es- 
timates of deformations in two dimensions and succeed in proving consis- 
tency within a subclass of deformations when observing random fields that 
are stationary but not isotropic. The subclass of deformations, however, is 
restrictive. In particular, it is not closed under post composition with rota- 
tions, which removes, in some sense, some complications for estimating the 
Jacobian generated by general deformations that can locally twist as well 
as stretch. On the other hand, Anderes and Stein [5] consider a large, non- 
parametric class of deformations. However, their results are methodological 



DEFORMED RANDOM FIELDS 



3 



Fig. 1. Left: One realization of a deformed isotropic random field. Right: The recovered 
isotropic random field using the estimated deformation from Figure 3. 

in nature and do not provide proofs of consistency. This paper contributes 
to bridging the gap between these two papers by considering the same flex- 
ible class of nonparametric deformations as in [5] while proving consistency 
for the estimated deformation as in [20]. For further references on densely 
observed deformed random fields see [9, 17, 29, 32, 33] and [35]. 

In this paper, we use kernel smoothed directional quadratic variations to 
estimate local properties of /, which are used to construct an estimate / 
of /. We establish sufficient conditions on the rate of bandwidth decay, in 
relation to the grid spacing, for strong uniform convergence of the kernel 
smoothed quadratic variations. There is a significant amount of literature 
studying the convergence of quadratic variations (see, e.g., [1, 6, 7, 8, 11, 15, 
18, 19, 24, 25, 30] and [41]). One of the crucial inequalities used in many of 
the recent convergence results is the Hanson and Wright bound for quadratic 
forms [21]. Indeed, we also depend heavily on this bound, and we use it to 
establish Claim 2 in Section 3.2, which, in turn, gives uniform convergence 
on compact subsets for the estimated complex dilatation and log-scale and, 
ultimately, the convergence of the estimated deformation. 

The kernel smoothed quadratic variations used in this paper are based on 
second order increments of the deformed process. Second order increments, 
rather than first order, are used in equation (4) to obtain sufficient spa- 
tial decorrelation for uniform convergence. Using higher order increments 
for quadratic variations is not new. They have been used in [24] and [7] for 
identification of a local fractional index and in [12] to identify the singu- 
larity function of a fractional process. The heuristic is that by increasing 
the order of increment, one can increase the rate of decay of the variance 
of the quadratic variation. However, this rate improvement holds only to 
a point, after which additional increments no longer improve the situation. 
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The interaction between the number of increments, the fractional index of 
the random field and the dimension of the domain of the random field is 
investigated in Chapter 3 of [3]. 

One of the main theoretical tools we use in this paper is the theory of 
quasiconformal maps. We believe this paper demonstrates how quasicon- 
formal theory can provide a flexible theoretical framework for estimating 
smooth invertible transformations, whereby making these objects available 
to statisticians for modelling a diverse range of physical phenomena. In two 
dimensions, an important object in the theory of quasiconformal maps is 
the complex dilatation : — > B (here is the observation region for Z o 
and D is the unit disk in the complex plane). A more detailed discussion 
of the complex dilatation is presented in Section 4. Besides characterizing 
quasiconformal maps up to post composition with conformal maps, the com- 
plex dilatation /i has two other useful properties. First, n can be interpreted 
as measuring the ellipticity and inclination of the local ellipse, which gets 
mapped to a local circle under the quasiconformal map that it character- 
izes. This is important to us for developing estimates of fi locally from Z o f. 
Second, the only requirement on is measur ability and ||/u||oo < 1- In other 
words, it suffices to measurably assign eccentricity and inclinations of local 
infinitesimal ellipses, and, by keeping the eccentricity bounded, there is a 
quasiconformal map that sends these infinitesimal ellipses to circles (unique 
up to post composition with conformal maps). This property allows us to 
find a smooth invertible transformation that corresponds to the estimated 
complex dilatation. 

The other object we use for estimating / is the log-scale r := log \ df\. We 
will discuss both /x and r in detail later. However, it is worth while to notice 
that /i and r provide enough information to uniquely specify the map / up 
to an overall rotation and translation. As we will see, one of the difficulties 
with r, as compared to /i, is that r lies in a complicated subspace of functions 
mapping M? to M. This is where we employ the Bergman projection as a tool 
to overcome the restrictive nature of the log-scale parameter r. 

Figure 2 illustrates the estimates fi and \df\ from the simulation shown in 
Figure 1. These are obtained by convolving the second order quadratic incre- 
ments with a Gaussian smoothing kernel and transforming these smoothed 
increments as discussed in Section 4. An estimated deformation / corre- 
sponding to fi and \df\ is shown in Figure 3 (left) along with the true defor- 
mation (right). The image shown in Figure 1 (right) shows Zo/o/~^, which 
"unwinds" the deformed process in an attempt at recovering the isotropic 
process Z. Note that the deformation / in Figure 3 is constructed using 
methods from both this paper and from [5]. To be explicit, all of the esti- 
mation methods from this paper are used for / with the exception of the 
Bergman projection, where the methods from [5] were used. The reason 
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Fig. 2. The estimated dilatation fi (left) and scale \df\ (right) using kernel smoothed 
second order directional increments and the results of Section 4- 



Fig. 3. Left: The estimated deformation recovered from the estimated dilatation and 
scale shown in Figure 2. The deformation was constructed from the estimated dilatation 
and scale using methods outlined in [5J. Right: The true deformation f . 



for this is that the computational techniques are not yet developed, to the 
authors' knowledge, for accurate approximation of the Bergman projection 
used in constructing /. However, since the Bergman space is a reproduc- 
ing kernel Hilbert space, there is potential for accurate approximation using 
spline methodology. 

There are three main parts to this paper. Section 2 discusses the as- 
sumptions on Z and the smooth invertible map /. In Section 3, we show 
that kernel smoothed directional quadratic variations converge uniformly on 
compact subsets with probability one. These directional quadratic variations 
are then used, in Section 4, to get estimates fi and f := log \ df\. Finally, in 
Section 5, we show how to convert ft and f to an estimated map / on simply 
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connected subsets U such that U CU, and we show that / converges to / 
uniformly on compact subsets of U with probability one. 

2. The random field Z and the map /. In this section, we list our as- 
sumptions on the isotropic random field Z and the smooth invertible trans- 
formation /. This section starts with a brief discussion of our assumptions 
on the autocovariance function of Z. Then, a detailed discussion follows on 
our assumptions for the smooth invertible transformation / : M'^ ^ . 

We require the following three conditions on the isotropic random field 

Z: 

Rl. Z is a constant mean Gaussian process on with autocovariance 

i?(|t-s|) = cov(Z(t),Z(s)); 
R2. R{\t\) = R{0) - |t|"-Fo(|t|"+^), as |t| ^0 for some < a < 2, 7 > 0; 
R3. R is away from the origin and there exists a c > such that | R^"^^ (t) \ < 
£qj, g^Yl sufficiently small t> 0. 

The assumption R2 establishes the local quadratic variation behavior of the 
process Z to be similar to that of a fractional Brownian sheet with Hurst 
index a/2. Informally, the assumption R3 ensures that the second order 
increments of Z have spatial decorrelation like that of a fractional Brownian 
sheet. 

Remark 1. Most of the following results can be extended, with some 
additional technical assumptions, to a larger class of autocovariance func- 
tions by replacing the principle term \t\" in R2 with L(|t|)|i|", where L is 
a slowly varying function at 0. The main difference is that the quadratic 
variation process defined in equation (4) below will need to be normalized 
by n~"L(l/n) instead of n~". 

Remark 2. The class of autocovariances satisfying R2-R3 encompasses 
a broad range of random fields that are continuous but not differentiable. Ex- 
amples include the Matern autocovariance function with smoothness param- 
eter < < 1 (see [31, 39] and [42]) and the exponential family exp(— c|t|"), 
where a £ (0,2). One way to extend our results to random fields with 
higher order differentiability is to use quadratic variations of higher or- 
der increments of the deformed random field to obtain sufficient spatial 
de-correlation. In the interest of space, we only prove results for the non 
differentiable case. 

Our basic assumption on the smooth invertible map / is that there exists 
a local affine approximation. In particular, 

(1) /(t + h) = /(t) + j)h + o(|h|). 
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where Jj := (^(t))jj is the Jacobian of the map / at t. This local lin- 
ear behavior is important, since we do not have replicates of the deformed 
random field Zo and, therefore, most of the statistical information is con- 
tained in the local variation of the process Z o f. When / behaves locally 
like the Jacobian matrix transformation, the distribution of the random field 
Z o f(t + h), as h varies in a small neighborhood of the origin for a fixed 
t, behaves similar to that of Z{Jjh.). Therefore, one can hope to estimate 
parameters of the Jacobian Jj using the local quadratic variation of one 
realization of the process Z o f near t. Of course, higher order terms in a 
Taylor expansion may also be estimated; however, these presumably require 
smaller neighborhoods for accurate estimation. 

In addition to the local affine behavior of /, we will need extra smoothness 
conditions. We discuss the following three notions of differentiability for a 
map f '.U ^ V between planer open subsets U, V: Frechet, Gateaux and 
C^{U). For a point ^U, f is said to be Gateaux differentiable at xq in 
the direction h if the limit 

l-^ /(xo + gh)-/(xo) 

exists. A stronger notion of differentiability is Frechet differentiable. The 
map / is said to be Frechet differentiable at xq if there exits a continuous 
linear map r : ]R2 ^ M2 such that 

/(xo + h)-/(xo) = T(h) + o(|h|). 

If such a T exists, and / has continuous partial derivatives, then T is the map 
corresponding to left multiplication by the Jacobian matrix Jj " . Clearly, if 
a function / is Frechet differentiable at xq, then it is Gateaux differentiable 
at xq and the Gateaux derivative in the direction h is equal to <^j"h. 

The third notion of differentiability C^{U) is satisfied if the partials 
exist and are continuous on U. This notion of differentiability is different, 
in that it is not defined pointwise. The reason for this is that there is not 
much one can say about the local behavior near xq of a function where the 
partials exist. In fact, it may be neither Frechet nor Gateaux differentiable. 
If, however, we require the partials ^ to exits and be continuous on U, 
then this is enough to imply Frechet differentiability at all points in U. 
Even more, it is true that / is C^{U), if and only if x i— > is continuous as 

a mapping from U into the space of continuous linear functions on M'^ (see 
Claim 8.9.1 from [14]). 

Define the class of C^{U) diffeomorphisms to be the set of all continuous 
invertible maps /:[/— >M^, such that / is C^{U) and is C^(V), where 
V = f{U) is the range of /. By the Inverse Function Theorem, necessary 
and sufficient conditions are that / be invertible, C^{U) and det Jj ^ 0. 
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We write as short for C^(]R2). Now, every diffeomorphism is Frechet 
differ entiable, and so there is a Jacobian matrix Jj such that 

/(t + h) = /(t) + j)h + o(|h|). 

Moreover, the directional derivative in the direction 6, denoted dof, is J^uq, 
where uq = (cos 0, sin 0). 

In the following paper, we will restrict the definition of diffeomor- 
phisms to have det Jj > on , which characterizes the diffeomorphism to 
be orientation preserving. In some sense, this is a trivial restriction, since, 
when / is a diffeomorphism, either det Jj > everywhere or det J/ < 
everywhere. These are referred to as orientation-preserving and orientation- 
reversing, respectively (see page 10 of [28]). Finally, define C"'(C/) diffeo- 
morphisms to be the C^{U) diffeomorphisms with order r continuous mixed 
partials. 

We now list some consequences of a diffeomorphic assumption on /, 
which we use in the following proofs: 

Dl. / is a quasiconformal map on bounded simply connected domains (see 
the Appendix of [4] for a definition of quasiconformal maps); 

D2. suptge I '^^^"'"'^ — — > as e — > for every compact set Q; 

D3. For any vector h 7^ and compact set 0, there exists a constant c such 
that |a{,^'^^i2(|/(s) -/(t)|)| < c|s-t|"-^ for ah s,t G 6 such that s/t 

(note that R and a are defined in assumptions R1-R3, and d^'"^^ is 
defined in the next section); 
D4. For every compact subset 0, there exists constants ci,C2 > such that 
ci|h| < I Jjh| < C2|h| for all h and all x G 0; 

D5. I J*h|" is Holder continuous in t G G for any h and compact set 0. 

Note that D3, D4 and D5 are the only statements that depend on the 
extra assumption rather than the C^. The proofs of D1-D4 for dif- 
feomorphisms are included in the Appendix of [4]. Notice that D5 follows 
from D4 and the fact that Jj has components. 

3. Kernel smoothed squared increments. In this section, we study the 
convergence of the kernel smoothed squared increments of the deformed 
process ^ (x) := Z o /(x) observed on some dense grid in a bounded, simply 
connected open subset C M^. The asymptotic regime we consider is as the 
grid spacing goes to zero and the region 0, stays fixed, which is sometimes 
called infill asymptotics. 

For a fixed nonzero vector h G M^, let 

Ah^t) :=y(t + h)-y(t), 

A£^y(t) :=AhA--^y(t). 
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If t is near dil., computing AJ^ may require observing Y outside of 0,. There- 
fore, we will suppose that we observe Y onQ plus some points within a small 
distance from the boundary dil. Now, for a function of two variables F{s, t), 
let A[7'"''i^(s,t) := Ag'AgF(s,t), where A{^ acts on the variable s and Ag 
acts on the variable t. Define d\i := h ■ V to be the directional derivative 
in the direction h and d^'"'^ F{s,t) := 9j^^9jJF(s, t), where acts on the 
variable s and acts on t. The following notation will be used throughout 
this paper: 

(2) 5(t):=(8-2-+i)|j)hr, 

(3) 0„:=On{ZVn}. 

That is, On is the grid of spacing 1/n in 0,. 

Here is a summary of the results of this section. First, we show Lemma 1, 
which establishes that E(A^, y(t))^ ^ n~'^g{t). Motivated by this lemma. 



we estimate g{t) by locally averaging the squared increments (A^^^y(w)) / 



n 



for the points w near t then show that there is enough spatial decorrelation 
for convergence. To this end, define B^^bi^) : ^ M as 

1 „ /w-t\(AS,„>'(w))' 



Here, is a convolution kernel satisfying certain conditions stated below. 
This section then culminates with Theorems 1 and 2 concerning the uniform 
convergence of Bn^b and d^Bnfi for some u 7^ 0. In particular. Theorem 1 
shows that, under appropriate conditions, 

sup|Sn^b(t) -g{t)\ — >0 w.p. 1 
tee 

as n ^ 00, 6 — > and n~^b~^ for all compact sets C fi. Theorem 2 
shows that, with some additional smoothness assumptions and — > 0, 
the directional derivatives d^Bn^b converge uniformly on compact subsets 
w.p. 1 to dug- 

3.1. Assumptions on K and Vl. The assumptions on K are as follows: 

Kl. K has bounded, continuous first and second order mixed partial deriva- 
tives; 

K2. //A'(x)dx= 1 and //|x|A'(x)dx<cx). 

Notice that these assumptions imply that K is Riemann integrable, K and 
the first partials of K are Holder continuous. Finally, we assume: 

01. J7 is a bounded simply connected domain of M^. 

Notice that this assumption ensures that the number of points in 0„ is of 
order n^. 
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3.2. Strong convergence of Bn^b o^nd d^Bn.b- For the remainder of this 
section, let Y := Z o f , B^^b be defined as in (4), g be defined as in (2), 
On be defined as in (3) and set Xt := Bn^bi^) ~ E-Bn,6(t). In the fohowing, 
'universal constant' means any constant that does not depend on n, b, Q or 
the process (Xt)tGe- 

Lemma 1. Suppose R1-R2, 01 and f is a C'^ diffeomorphism. Then, 

E(A,\y(t))2 



sup 
ten 



9{t) 







as £ loo. 



Proof. By assumption R2 we can write R{\t\) = R{0) - + r{\t\), 
where r(|t|) = odtl^+T) as \t\ 0. Write E{Al^^Y (to))^ for to G H as a sum 
of two terms 

E(ALy(to))' = Af/){cov(y(s),y(t))}|,,t=to =Ii+l2, 

where 

(5) Ji = Af/){-|/(s)-/(t)r}|,,t=to, 

(6) T2 = Af/){r(|/(s)-/(t)|)}|,^,^,^. 

Write the increment operator A^j^ as a linear filter, so that its action on 
a function Q -.M? — > M can be expressed as A'^^Q{t) = J2'j=odjQ{t + esj), 

where dj = (— and Sj = jh. Now, the first term can be computed as 

2 

(7) j,/e" = _ ^ did,\ f (to + esi) - f (to + £s,r/e'' 

i,j=0 
2 

(8) =- ^ did,|jJ°(s,-s,)r + o(l), 

i,j=0 

where o(l) — > uniformly over to G 17 as e — > by D2 and D4. 
Similarly, the second term can be computed as 

(9) T,/." = f ^^rf /(l/(to + ^^^);/(to + ^^.)l) ^ 

i,j=0 

where sup^^^^^ |r(|/(to + esj) — / (to + esj)|)/e"| converges to zero by R2, D2 
and D4. Combining terms Ii and X2, we get 



i?(A,h>^(to))V^° ^ - E d^dj\J}"{s, - s,)r 



i,j=0 
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uniformly for to € ^2. This completes the proof, since — J2jj=o didj\J^° {si — 
SjT = g{to). □ 

Lemma 2. Suppose R1-R3, 01 and f is a diffeomorphism. Then, 
there exists a constant c > such that 

|EA2/„y(t)A2/„y(s)|<cn-4|^-sr-4 

for all s,t £ $7 such that |s — t| > |3h/n|. 

Proof. The idea is that 

|EA2/„y(t)A2/„y(s)| = |Ag>(|/(s) - /(t)|)| 

«|n-^4''')ii(|/(s)-/(t)|)| 
< cn |s — t| . 

To make this precise, let F(s,t) := cov(y(s), y(t)) = R{\f{s) - /(t)|) and 
H be the 2 by 2 matrix with each column h/n. Then, 

|EA2/„y(t)A^/„y(s)| = |Ag'jF(s,t)| 

_4 /• f d^^''^^F{s + Htt + Hri)d^dri 



n 



'[o,ip J[0,1] 

<cin~^ I I \s-t + H(^-ri)r-^d$dT] 

i[0,l]2 J[0,l]2 



< C2n ^ sup |s - t + h(^i +^2 - ??i - '?/2)/'rar ^ 



= C2n ^ sup |s — t + 2hr/n|" ^ 

-l<r<l 

< C3n~^|s — tl"^"^ when |s — t| > |3h/n|. 

Notice that the above proof requires that \d^''^^ F{s,t)\ < ci\s — t|°~^ 
which is why we need D3. □ 

Claim 1. Suppose R1-R3, K1-K2, 01 and f is a C'^ diffeomorphism. 
Then, 

sup|E(B,,b(t))-<7(t)H0 
tee 

for all compact subsets Q CQ as oo, 6— > and n~^b~^ = o(l). 
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Proof. First, let 5'n(w) := n°E(A^y^y(w))^ for w G il. We show 



E 



w — t ^ 
K[ —— ]g{w)+ei 



= g{t) + em + en + ei, 

where ei = o(l), en = 0(n^^&~'^) and em = o(l) uniformly on compact sub- 
sets of O as n ^ oo, 6 ^ and n^^b^^ 0. 

To show the results about ei and en, we need to control the error when 
approximating Riemann integrals of Holder continuous functions on by 
Riemann sums on This error is bounded by the difference between the 
upper and lower Riemann sums, which is bounded by a^an~^, where a is 
the Holder constant of the function and is a constant only depending on 
the region $7. To show ej = o(l) and en = 0{n~^b~^) uniformly for t G 0, 
we will use the fact that b~^K{{- — t)/b) is Holder continuous with Holder 
constant cb~^ for some constant c. 

To show ei = o(l), fix some compact subset C and notice that 

Iei|<i2 E \K\(^)\9nM - g{^)\ 



n^b"^ „ 



< sup |5„(w) 




1 




The term sup^^^^ |<7n(w) — g{'w)\ — > by Lemma 1. Now, by the comments in 
the previous paragraph, the Riemann sum ^753- X]wen„ l-^l(^T^) approx- 
imately / |iir|((x — t)/6) dx (which is bounded) with error 0{n~^b~^). 
Therefore, ej ^ as n — > 00, 6 — > and n~^b~^ 0. 

Similarly, to show en = 0{n~^b~^), we notice that the Holder continuity 
of K and g are sufficient for the Riemann sums of b~'^K{{- — t)/b)g{-) to 
converge to the Riemann integral with an error en = 0{n~^b~^) uniformly 
in t e e. 

Finally, to show eni = o(l), we need that 
1 

62 



K[—-]g{i^)di^= / K{^)g{b^ + t)d^^g{t) 

n \ b J J J{Q^t)/b 



as 6 ^ uniformly in t G 0. Here, we use the Holder continuity of g and 
assumption K2. Notice that the error term eni does not converge to zero 
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uniformly in t € 17. This is why we can only show the result uniformly on 
compacts instead of uniformly on $7. □ 

In what follows we will not only show convergence results about Bnfi{t) 
but also 9u-Bn,b(t), and Bnfi{t) — Bn,b{^) all of which have the form 

(10) Q„ = n""2^ (A2/„y(i))2/Ci, 

ien„ 

where 1C\ may depend on t,s G 17 and the bandwidth parameter h. 

Claim 2. Let Qn be defined as in (10). Suppose R1-R3, 01 and f 
is a C"^ diffeomorphism. In addition, suppose that there exists a function 
G(t,s,6), such that |/Ci| < G(t,s, 6) for all s,t G $7 and b in a neighborhood 
of the origin. Then, for all e > 0, b sufficiently small and n sufficiently large, 

r^r,^ 1 / C2en^ c^e^n^ \ 

P[|<3„ - Ee,.| >el < ci»p(-^j|^ A ^ji^), 

where ci, 02,03 are universal constants. 

Proof. Writing /Ci = IC^ — /Cr, we get the decomposition Qn = Qn — 
Qn, where 

Q^=n°-2^(A2/„y(i))2/C+, 
Q^=n-2;^(A^/„y(i))2/c^. 

Therefore, 

P(|Qn - EQn\ >e)< p(\Qi - EQi\ > |) + P (\Ql " ^Ql\ > • 

First, we find a bound for P(|Q^ — EQn\ ^ §)• Let Ay be the column vector 
with elements n"/2-i^2^^y(i)^+/2 f^^. j ^ ^^^^^ ^+/2 _ (/c+)i/2. Let 

S = EAYAy^ be the covariance matrix for AY, so that Qj^ = AY'^AY = 
W'^T.W, where W is a vector of i.i.d. standard Gaussian random variables. 
Also, let S(i,j) denote the matrix entries of T, for i,j E 

Using the bound on quadratic forms for Gaussian random variables found 
in Hanson and Wright [21], we now get 

P[\W^J:W -EW^j:W\>e]<2exp(--^^^ A-^^^), 

\ ll^abslb ||Sabs|li?/ 
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where || • II2 and || • \\f are the spectral and Frobenius matrix norms, respec- 
tively, and Sabs is the matrix with elements |S(i,j)| for i,j G 0„. Now, 

^ 1 1 — 4i- "la— 4 

< C3 TT-^ — n 1-1 

for all |i — j| > |3h/n| by Lemma 2. By assumption, \K^^'^K^^'^\ < G(t,s,5). 
Therefore, for all i, j G such that |i — j| > |3h/n|, 

(11) |S(i,j)|<C3G(t,s,6K-6|i-jr-4. 

To finish the proof, we show ||5]abs|lF ~ 0(|G'(t, s, 6)pn~^) and ||Sabs||2 = 
0(|G(t, s, 6)|n~^) uniformly for all t,s E 0, h sufficiently small and n suffi- 
ciently large. 

To show the bound for ||Sabs|||'5 notice that 

Es(i,j)2< S(i,j)VciGV""« 

ijen„ |i-j|<|3h/n| |i-j|>|3h/n| 

= * S(i,j)2 + ciG2n2"-80(n-2a+6) 

|i-j|<|3h/n| 

where the last equality is because < maxiS(i,i) = 0(n~^G(t, s, 5)) 

by Lemma 1. To get =*, there are some technical difficulties, but the heuris- 
tic is, when < a < 2, 

Y n-^|i-j|2"-«xc4 / Ixp^-^dx 

|i-j|>|3h/n| ^{l/n<|x|<l} 

= C5 t'^^-Ut 
Jl/n 



0{n 



-2a+6\ 



For the full details, see [3], Lemma 3, page 41. Finally, to show the bound 
for II Sabs II 2) notice that 

||Sabs||2<max^ |S(i,j)| 



ii-jI"-' 



<max Y |S(i,j)|+C3Gmax Y „6-a ^y (11) 

^ "|i-j|<|3h/n| '*^''"|i-j|>|3h/n| " 
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= 0(n-2G') + C3Gn°-Vax ^ n-2|i-j|"-^ 

|i-j|>|3h/n| 

= 0{n-'^G)+0{n-^G), 

where the last equahty uses the fact that < a < 2. 

This estabhshes the desired bound for P{\Qli — ^Qn\ > §)• The result for 



P(IQn ~ ^Qn\ ^ |) is exactly similar. This completes the proof. □ 



Corollary 1. Fix a point tg G i^, suppose R1-R3, 01 and let f be a 
diffeomorphism. If K is bounded, then 

\Bn,b{'to) — E-Bn,f,(to)| -^0 

as CO, 6 — > and n~^b~'^ = 0{n~^) for some P > 0. 

Proof. This follows by Claim 2 and Borel-Cantelh using /Ci = ^Er(i^) 
and G(t,s,6) = 6-2||is:||^. □ 

Corollary 2. Fix a point to G ^, suppose R1-R3, 01 and let f be a 
diffeomorphism. If K has continuous partial derivatives and u / 0, then 

|9uS„,fe(to)-EauS„,b(to)|^0 

as n—>- oo, b—>^0 and n^^b^^ = 0{n^^) for some /3 > 0. 



Proof. This follows by Claim 2 and Borel-Cantelh using /Ci = ^ x 

b 



{duK){^) and Git,s,b) = b-^duK\\^. □ 



The following corollary will be used for the uniform convergence of to g 
in the next subsection. Remember that Xt is defined as Bn.b{^) — Ei3„^f,(t). 

Corollary 3. Fix a point to G f^, suppose R1-R3, 01 and let f be a 
diffeomorphism. If K is Holder continuous, then 

(12) P(|Xt -Xsl >e) <ciexp[^ — |^ _ ^| A |^_g|2 

where ci,C2,C3 are universal constants. 

Proof. First, write Xt — Xs in the form Qn — ^Qn, where Qn '■= -Bn,fe(t) ■ 
Bn^bi^)- Then, the corollary follows by Claim 2 using 

so that |/Ci| < c6~^|s — 1| for some Holder constant c > 0. □ 
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3.3. Uniform convergence of Bnfi and d^Bn^b- In this subsection, we use 
the results from the previous section to estabUsh the uniform convergence 
of Bn^b and the directional derivative d^Bnfi on compact subsets of the 
observation region $7. These results are then used to establish consistent 
estimators of the complex dilatation fi and log-scale r of the diffeomorphism 
/ in Section 4. 

Lemma 3. For any b> and a>e, we have 



(ae""* Al)dt<2J^- and / (ae""* A 1) < 2 



V Jo 



b 



Proof. Given 6 > and a > e, let 7 = \/b^^ logo. Then, 
r {ae-'''' M)dt = j+ r ae-^*' dt 

Jo J"i 



<7 + / -ae"^* dt 

h 7 

< 27 (since a > e =^ 7 > I/267). 
Similarly, putting v = b~^ log a, 

/ {ae-^^ M) dt = V + / ae'^'Ht = v + - <2v. 
Jo Ju 

This completes the proof of the lemma. □ 

Lemma 4. Let Q be a compact subset of il.. Suppose the assumptions 
in Corollary 3. Let ti, . . . ,tm and si, . . . , be arbitrary points in 0, where 
m > 2. Let M = maxi<j<m l^ti — -^si I a'^d 5 = maxi<j<m |ti — Sj| . Then, for 
any r > 0, 

/ C2rn%^\ ( c-^r'^TV^b^ 

P(M > r) < c\m exp I 1 + c\m exp 



where c\, C2 and C3 are universal constants. 

Proof. By Corollary 3, for each r > and t,s G 9, 
, N ^ corn'^b^ c^r'^nH 



|t-s| |t-s|2 

<ciexp(-^^j+c,expf---— 
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From the above bound, we see that 

P[M >r)< cimex.pl j + cimexpl ^ j . 

This completes the proof. □ 

Claim 3. Fix a compact set Q C il, and a point to £ @. Let M = 
maxtee l^t ~-^tol ^'^'^ suppose the assumptions in Corollary 3. Then, there 
exists universal constants Li and L2 such that, for all R> 0, we have 

p(m> <Liexp{-L2mm{Rn,R^}). 



Proof. Suppose we have a sequence of finite sets Aq, Ai, A2, . . . C Q 
and constants c<l < B and D satisfying the following properties: 

(i) Ao = {to}; 

(ii) Ak C Ak^i for all k; 

(iii) l^fcl < for all k, where |^;;| denotes the cardinality of Af^; 

(iv) For each k>l and each t £ A^ there exists a 'parent' tp G Ak-i such 
that |t — tpl < Dc^ (note that t = tp if x G Af^_i); 

(v) The sequence has a 'limiting denseness property' in the sense that, for 
any nonnegative continuous function / on Q, we have maxtge/(t) = 
limfc^oo maxteAfe /(t). 

It is easy to see how the constants c, B and D can be chosen, with D oc 
diam(0), and the sets {Afc}fc>i constructed by successive dyadic partition- 
ing. For each A; > 1, let 

Mfc = max|Xt -^t„|- 

teAfc " 

Applying the 'limiting denseness property,' we see that M < J2'h=i -^fc- Foi" 
each /c > 1, let = sr^c^ . j . Then, rk > and J2k'^k = 1- Thus, for any 

R>0, 

P (m > < P f Mfc > ^ for some A; > 1 

nb'^ I V nh^ 



(13) 

By the tail bound in Lemma 4, we have 



< ci-B''(exp(-C5i?/cn) + exp(-C6-R A:^)) 
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where C5 and cg are universal constants. Thus, if ii > cjlogB for some suit- 
ably large constant C7 (that need not depend on n), then, from the above 
bound and (13), we can conclude that 

p(m> -^t^ < C8 exp(— cgi^n) + C8 exp(— cioi?^). 

The condition R > cj log B can be dropped by choosing large enough, 
because B, cg and cio do not vary with n. This completes the proof. □ 

Now, we obtain Theorems 1 and 2, which are the main results of this 
section. These establish uniform convergence results for the kernel smoothed 
squared increments. 

Theorem 1. Suppose R1-R3, K1-K2, 01 hold and f is a C'^ diffeo- 
morphism. Then, for all compact sets Q Cfl, 

sup\Bn,b{t) - g{t)\ — >0 w.p. 1 
tee 

as n ^ 00, 6 — > and n~^b~'^ = 0{n~^) for some (3 > 0, where g(t) := (8 — 
2"+i)|J}h|°. 

Proof. Let be a compact subset of Q and let to G ©. Let 6„ — > as 
n ^ 00, so that n~^b~^ < cn~^. Set M„ = sup^ge l^t — -^tol) where Xt := 
Bn,bAt)-EBn,b.M. Then, 

SUp|S„,;,Jt) - (7(t)| <SUp|5„,;,Jt) - E5„,;, Jt) | + SUp | EB„,b„ (t) - ^(t) | 

tee tee tee 

< Mn + \Bn,b„{to) - ES„,b„(to)| + 0(1) by Claim 1. 

Now, by the bound in Claim 3, we have 
00 00 

PiMn > < L2exp{-L3mm{n^+^/^,nl^}) < 00. 

n=l n=l 

Therefore, M„ — > with probability one. Finally, n~^b~^ = 0{n~^) for some 
/? > is a sufficient condition for |Sn,6„(to) — Ei?„ ;,,^(to)| by Corollary 
1. □ 

Theorem 2. Suppose the assumptions of Theorem 1 hold along with the 
additional assumption that K is a compactly supported kernel and f is a 
diffeomorphism. Then, for all compact sets G C O, 

sup |<9u-B„,6(t) - dug{t)\ — > w.p. 1 
tee 

as n ^ 00, 6^0, n~^b~^ = 0{n^^) for some /J > such that (3 > 1 — 47. 
Note that ^ > is the number appearing in assumption R2 and g{t) := 
(8 _ 2^+1)1 jth|«. 
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Proof. Let be a compact subset of Q and let tg G 0. Let 6„ — > as 
n ^ cxo so that n~^b~'^ < cn~^ . Set M„ = sup^ge l^t — -^tol) where Xt := 
9uSn,b„(t) - E5u5„,fe„(t). Then, 

sup|auS„,b„(t) - 9u5'(t)l < sup|5u-B„_fe„(t) - E(9u-B„.6„(t)| 
tee tee 

+ SUp|EauBn,6„(t) -5u5'(t)| 

tee 

< M„ + |auS„,;,Jto) - EauB„,f,„(to)| +0(1), 

where o(l) is estabhshed by similar methods to that of Claim 1 using the 
fact that ii' is a compactly supported kernel. 

Now, Corollary 2 establishes that |c?h-Bn,6„(to) — Ec)h-B„^;,^(to)| 0. The 
proof that Mn — — > is similar to Theorem 1, with the exception that, now, 
instead of the bound from Claim 3, one can show 

P(^M„>-^^ <Liexp(-L2min{i?n,i?2})^ 
which is sufficient to prove the claim. □ 

4. Consistent estimates of (/x, t). At this point, it becomes advanta- 
geous to switch to complex variable notation, so that points in the plane 
(x, y) S correspond to points in the complex plane a; + £ C. Under this 
correspondence, diffeomorphisms f{x,y) = {u{x,y),v{x,y)) of can be 
considered as diffeomorphisms of C by writing f{x + iy) = u{x,y) + 
iv{x,y). For the remainder of the paper, we use Re / and Im/ to denote the 
real and imaginary parts of the complex representation of the map /. 

The main utility of switching to complex notation is that we can directly 
use the results and techniques of conformal and quasiconformal theory to es- 
tablish consistent estimates of /. This section starts by defining the complex 
dilatation fi and log-scale r of a diffeomorphism. We then conclude by 
showing how the results of the previous section can be used to construct con- 
sistent estimates of fi and r. In the next section, we show how the estimates 
of /.i and r can be used to establish consistent estimates of /. 

For a function f £ C^{U), define the complex derivatives 

2\dx dy )' 2\dx dy 

The complex dilatation and the log-scale are then defined as 

(14) f^-=df/df, 

(15) r:=log|a/|. 

The complex dilatation ^ characterizes the infinitesimal ellipse with inclina- 
tion arg(— /i/2) and eccentricity that gets mapped to an infinitesimal 
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circle under the image of /. In addition, ^ uniquely determines / up to 
post composition with conformal maps. The log-scale r is then used to re- 
cover the conformal post composition, so that, together, /u and r uniquely 
determine / up to a rotation and translation. For a short introduction to 
quasiconformal theory, see the appendices of [5] or [4]. For more a complete 
treatment, see [2, 26, 27] and [28]. 

In the previous section, we constructed a sequence of functions i?„^ft^(t), 
which converge uniformly on compacts as n — > oo to (8 — 2"^^)| J*h|", where 
h := (/ii,/i2) is a vector of our choice. Using complex variable notation, we 
can write | J*h| = \hdf + hdf\, where h = hi + ih2- By factoring out 9/, 
we get \hdf + hdf\°^ = \df\°'\h + h(j\°^ . Now, by choosing h = + i (for 
increments in the north-south, east-west and diagonal directions), we can 
construct three functions Wi^n, W^2,n and W^^n, which converge to 
— ;u| and |f?/||l + i — respectively. In particular. 



W^l,n(t) 



nil 



4(? 



W^3,n(t) 



nil 



where the convergence is uniform in t on compacts subsets of as n — > oo. 
Notice that Wi^n, W^2,n and W^^n are the factors of stretching that the affine 
transformation Jj applies to the lines in the horizontal, vertical and diag- 
onal directions. Therefore, the points, (VF{"^,0), (0,VF2~^) and (1^3~^, ^^3"^) 
asymptotically lie on the ellipse that gets mapped to a circle with unit radius 
under the affine transformation induced by the matrix Jj . Since the general 
equation for an ellipse is ax^ -|- bxy + cy'^ = 1 , we have 

On := Wi n^a, 

Cn ■■= c, 

K:=Wl,,-Wl^-Wl^^b. 

The area of the ellipse specified by ax"^ + bxy + cy"^ = 1 is 2Tr/\/4ac — 6^. 
Since Jj sends the ellipse ax^ -|- bxy + cy"^ = 1 to the unit circle, we have that 
TT = det( J)) ^ , which gives 



AanCn -bl^ 2det(j)) = 2{\df\'' - 

Now, to solve for (/i, r), first notice that -|-Cn — > -|-/Up -I- |1 — /i| 

2|a/|2 + 2|5/|2. Therefore, ^AanCn-hi + an + Cn^4\df\\ Similarly, 

(16) 4Re(;u) = |l-h/ip-|l-//p, 

(17) 4Im(;U) = |l + i + /i(l-i)|2-|l + ^|2-|l-;^|2^ 
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which gives 

(18) ^ ^Re(^), 

(19) % ^ ^Mm), 



(20) log(^4a„c„ - 62 + a„ + c^) — ^ 2t + log 4. 

Therefore, under the conditions of Theorem 1, we can construct estimates 
fL and f that converge to // and r, respectively, where the convergence is 
uniform on compact subsets of with probability one as n — > cx). Moreover, 
under the extra conditions of Theorem 2, dfi dfj, uniformly on compact 
subsets of with probability one. 

5. Estimating /. In this section, we show how to construct an estimate 
/ on a simply connected domain U such that [/ C O. Then, we show that / 
converges to / uniformly on compact subsets of U. The construction of / is 
on U instead of because we need the uniform convergence of fi, dfi and f to 
establish the convergence of /. It is open as to whether one can construct / 
on the full observation region that converges uniformly on compact subsets. 

We start, in Section 5.1, with a discussion on how / can be recovered, 
uniquely up to a rotation and translation, from the true fi and r. This will 
indicate how we recover / from the estimated fi and f. Finally, we show 
f ^ f uniformly on compact subsets of U in Section 5.3. 

5.1. Recovering f from (fi, r) . First, let U he a simply connected domain 
such that U cQ. The diffeomorphism / now satisfies f = g o on U , 
where is the unique normalized quasiconformal map with dilatation /i 
that maps U to the unit disk D (see the Appendix in [4]). Since / and 
have the same complex dilatation, g = f o is a conformal map defined on 
D. It turns out that this decomposition of / is useful, in that the complex 
dilatation n determines and r is used to recover the conformal map g. 

To see how to recover g from r, notice that df = d{g o /^) = (^g' o /^) 5/^. 
Therefore, 

(21) log W\ = log o /-I - log 19/^1 o f-\ 

Since g is conformal on D, log (7' is holomorphic on B. Moreover, \ogg' = 
^og\g'\+iavg{g'). Therefore, using (21), Relogsr' = r o /"^ - log o 
which can be recovered from (/i,r). Since the real and imaginary parts of 
holomorphic maps are harmonic conjugates, which are unique up to a con- 
stant, we can recover Imlogf?' + 0, where ^ G M is an unknown factor. Now, 
by exponentiating, we can recover e'^^g' . Then, 



e^^g{z) + c= I e'^g'{w)dw 
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where the integral is taken over a hne connecting zq to z. Therefore, fi and 
T are sufficient to recover f = g o on U up to a rotation and translation. 

5.2. Constructing f from {fi,T). The technique for recovering / in Sec- 
tion 5.1 would work if we knew that there existed a quasiconformal map 
/ with complex dilatation fi such that f = log \df\. Unfortunately, there is 
no simple condition on f for the existence of such an /. The main problem 
is that we do not precisely measure log|(/'|, which is required to be har- 
monic. Instead, we only have an estimate of loglt^'j. The estimate, which is 
motivated by (21), is defined by 

(22) l^h=To/ri-log|a/^|o/ri. 

Since log \g'\ is not guaranteed to be harmonic, it may not always be possible 
to find the harmonic conjugate used to recover logg' . In what follows, we 
notice that log g' is in the Bergman space of holomorphic functions with finite 
L2 integrals. We then use the Bergman projection to find a holomorphic 
function whose real part approximates logj^'j. 
We define our estimate / of / in the region U as 



(23) f = goff 



A" 



where is the unique normalized quasiconformal map sending [/ to B, and 
the function g is the holomorphic map, unique up to translations, defined 
on the unit disk B and satisfying 

g' = exp{Vlog\g'\), 
where the operator V is defined by 

rh(w) :=- f , dxdy -Reh(0), 

vr (1 - zwy 

where z = x + iy. The integral transform in the above definition is the 
Bergman projection (see [16] for an introduction to Bergman spaces). 

To motivate our choice of operator V, we first mention that the true 
conformal map g satisfies 



[V\og\g \){w) = V\ ]{w) 



- [ dxdy + - [ dx dy - log \g' 

vr Je) (1 — zw)"^ vr Je) (1 — zw)"^ 

log g'{w)+ log g'{0)-log\g'{0)\ 
log g' (w) + i9, 
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where 9 = — Imlog (7'(0). In the above computation, we used the fact that 
for any conformal map g defined on D, the holomorphic function log 5' is in 
the Bergman space ^^(B) (this is true by Theorem 9.4 of [36] along with the 
fact that the Bloch space is a subset of the Bergman space). Therefore, the 
projection V can be used to recover the harmonic conjugate of log|(7'| up 

to an unknown constant 6. In what follows, we show log|(7'| ^logl^r'l and 

■piog Iff'! ^ "Plog \g'\ uniformly on compact subsets of ID as n ^ 00. 

5.3. / converges to f. We show that, under appropriate conditions, / 
converges uniformly on compact subsets of U, with probability 1. First, we 
establish the following lemma. 

Lemma 5. Suppose that fin,fJ' £ C'^{U) are complex dilatations such that 

IJ-n ^^^^ fi and d^n ^^^^ dfi on a bounded, simply connected domain U . 
Suppose, in addition, that one can extend iJ,n,fJ- to functions fJ-n^fJ'* G C^(VF) 

on a simply connected domain W containing U such that fi and 

d^n ^'^^'^ dfi. Then, 

(24) iog|a^jo/-i^iog|a^|o/;i, 

(25) V log \df,„ I o /-i ^ V log la/^i o f-\ 

uniformly on compact subsets o/D as n — > 00. 
Proof. First, decompose /^^^ and so that, 

ffin=hn°fn, ffi = hof, 

where /„ and / are normalized quasiconformal maps on the whole plane 
with complex dilatations obtained by extending /i„ and fx to the whole 
plane by smoothly truncating to zero away from U. Here, /i„ and h are 
conformal maps sending fn{U) and f{U) to the unit disk D, respectively. 
Note that the truncation must be done in such a way that jln has uniformly 

bounded compact support, is as smooth as jln fi and djln djl 
(existence guaranteed by the existence of the extensions , /i* on W) . Now, 

the quasiconformal maps /„ converge uniformly to / on U , since /i„ /2 
with uniformly bounded support (see Lemma 1 on page 55 of [2]). 
Notice = {h'n o /„) 5/„,, which gives 

(26) (log9/;.J o = logK, o /i^i + log5/„ o 4-1. 

See the Appendix in [4] for a discussion on how to define a continuous version 
of log9/^„ and log9/„. To establish (24) and (25) we show that both terms 
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in (26) converge uniformly on compact subsets of D, as well as the result of 
applying the operator V to both terms. 

For the first term, \ogh'^ o , in (26), notice that = fn° f^^ and 
h^^ = f o ff^^ ■ Since fn^f uniformly on [/, / is Holder continuous on U 
and f^^ ~^ ffj^^ uniformly on compact subsets of B, ^ h^^ uniformly 
on compact subsets of ID (this follows by Corollary 9 and Lemmas 10 and 
13 of [4]). Since the functions := and ^ := are conformal maps of 
D, log^^ and log^' are both holomorphic and log^^ logC uniformly on 
compact subsets. Noticing that log^^ = — log/i'„ ° and log^' = — log/i' o 
, gives 

(27) logh'^oh-^^logh'oh"^ 

uniformly on compact subsets of D as n — > oo. In addition, log^^ and log^' 
are both in the Bergman space ^2(15) (and are therefore unaffected by the 
Bergman projection), which establishes that 

(28) V\ogh'^oh~^ ^Vlogh' oh^^ 

uniformly on compact subsets of D as n ^ oo . 

For the second term log9/„ o in (26), notice that the results in 
the Appendix of [4] establish that log 5/ is Holder continuous on U and 

logd fn logdf. Therefore, logd fn o /~J converges to log 9/ o uni- 

formly on compact subsets. Moreover, since the continuity of log 5/ on C 
implies that it is bounded on U, logd fn ° fj^^ also convergences in L2{p). 
Therefore, 

(29) log dfn o ^ log df of-\ 

(30) r log dfn o ^ r log df o /-I 

uniformly on compact subsets of D as n — > oo. The last convergence is due 
to the fact that the Bergman projection is a bounded operator on L2(0) and 
that convergence in the Bergman space implies convergence on compacts. 
Finally, (27), (28), (29) and (30) estabhshes the lemma. □ 

Theorem 3. Suppose R1-R3; 01, K1-K2, K is a compactly supported 
kernel, f is a diffeomorphism, n ^ oo, 6^0 and n~^b~^ = 0{n~^) for 
some /3 > such that /3 > 1 — 47. Let U be a simply connected open subset 
of the observation region O such that U Cft. Then, the estimated map f, 
defined on U by {23), converges to e^^f + c uniformly on compact subsets of 
U with probability one, where 9 is an unidentifiable rotation angle and c is 
an unidentifiable translation. 
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Proof. The results of Theorems 1 and 2, along with the comments 
made in Section 4, establish that /i — > /U, dfi and f — > r uniformly on 

U with probability one as n — > oo. Therefore, — > uniformly on compact 
subsets of U (by Corollary 9 of [4]). It is now sufficient to show that g 
converges uniformly on compact subsets of B to g (where sufficiency is by 
Lemma 11 of [4]). 

We first show 

(31) log\g'\ ^log\g'\, 

(32) ri^i^riogig'i 

uniformly on compact subsets of B as n — > oo. Remember, log \g'\ is defined 
by 

(33) l^|:=ro/.i-log|a/^|o/.i. 

Lemma 5 immediately establishes the required convergence for the second 
term log |(?/^| o f7^. The first term f o f7^ converges to r o both uni- 
formly on compacts of B and in L2(B). This follows by Lemma 13 of 
[4], since f j that r is Holder continuous and bounded on U 

(since / is assumed to be a diffeomorphism) . Since the Bergman pro- 
jection is a bounded operator from L2(B) to yl2(B), we also have that 

Vt o fji^ 4^0^) -p^ o /~^. The convergence is also uniformly on compact sub- 
sets by Lemma 12 of [4]. This establishes (31) and (32). 

Now, by the comments made in Section 5.2, "Plog \g'\ = \ogg' + i6. There- 
fore, by (31) and (32), exp(7^1og l^'l) converges uniformly on compacts to 
e^^g'. Since U is simply connected, exp(Plog l^'l) has an antiderivative g 
such that g' = exp(Plog j^'l) and g — > e^^g + c uniformly on compact subsets 
of B. Therefore, / converges to e^^ f + c uniformly on compact subsets of U. 
□ 
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